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, Thermal and chh'al critical exponents of the fully frustrated XY model on 

^ , a square-lattice are obtained from a finite-size scaling analysis of the free 

O 

^ ' energy of chiral domain walls. Data were obtained by extensive Monte Carlo 

^ I transfer matrix computations for infinite strips of widths up to 14 lattice 

, spacings. Two transfer matrices were implemented, one for each of the two 

principal lattice directions. The results of both are consistent, but the critical 

exponents differ significantly from the pure Ising values. This is in agreement 

with other recent Monte Carlo simulations. Our results also support the 

identification of the critical behavior of this model with that along the line of 

transitions of simultaneous ordering or becoming critical of Ising and planar 

rotor degrees of freedom in the XY-Ising model studied recently. 
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I. INTRODUCTION 



There has been considerable interest, both experimentally and theoretically, in phase 
transitions of two-dimensional, fully frustrated XY models. Most studies have been moti- 
vated by their relevance to Josephson junction arrays in a magnetic field, where these models 
describe the superconducting-to-normal transition at half a flux quantum per plaquette0, but 
there are also interesting theoretical questions regarding the identification of the universality 



In the standard XY model without frustration the relevant symmetry is the continuous 
U{1) symmetry which, as is well-known, is responsible for a Kosterlitz-Thouless phase tran- 
sition. The fully frustrated XY model has a richer behavior with a low-temperature phase 
with critical fluctuations associated with the U{1) symmetry, accompanied by a broken, 
discrete Z2 symmetry. In principle, there are two different ordering scenarios: ordering can 
take place in two stages via separate XY and Ising transitions, or both symmetries can be 
broken or become critical simultaneously, which yields a single transition, presumably in a 
new universality class. 

On a square lattice the model can be defined by the Hamiltonian 



where = J (J > 0) for horizontal rows and Jij = ±J for alternating vertical columns. 
Owing to the presence of an odd number of antiferromagnetic bonds in each plaquette, the 
model is fully frustrated. This leads to a double degeneracy of the ground state, which is 
of course accompanied by an additional continuous degeneracy, a manifestation of the U{1) 
symmetry. One can introduce an Ising-like order parameter Xpj the local chirality, which 
measures the sense of rotation of a spin of unit length s = {cos 6, sinO) along the sides of a 
given plaquette p. In the ground state there is antiferromagnetic order of the local chiralities 



Early Monte Carlo simulations led to the conclusion that the chiral order parameter of 
the fully frustrated XY model has pure Ising model critical exponentsi'Si'iffli, but more 
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estimates of the exponents have differed significantly from the Ising values. 
In particular, predictions were made for the values of chiral critical exponents of fully frus- 
trated XY model. These predictions are based on results for the XY-Ising mo del0'0, which 
is expected to describe the critical behavior in these systems. These ideas are supported 
by recent Monte Carlo simulationJlllli. However, in view of results that suggest pure Ising 
critical exponents, and in the absence of precise agreement among the more recent esti- 
mates, the current state of affairs is unsatisfactory from a computational point of view. The 
additional numerical results presented in this paper may serve to help settle the issue. 

We report results of extensive numerical calculations of the chiral domain wall free energy 
of the fully frustrated XY model in an infinite strip geometry. Because of the continuous 
nature of the phase variables of this model, the transfer matrix is intractable for numer- 
ically exact computation of its eigenvalues. We therefore use the Monte Carlo transfer 
matrix methodEl to obtain the free energy from the largest eigenvalue of the transfer ma- 
trix. Exploiting the anisotropy of the hamiltonian (|1|) we use two different implementations, 
a "horizontal" {i.e. column-to-column) and "vertical" {i.e. row-to-row) transfer matrix. 
These approaches yield results in reasonable mutual agreement. From a finite-size scaling 
analysis of data for strips of widths up to 14 lattice spacings, we have estimated the thermal 
exponent u, the exponent rj associated with the Ising-like order-order correlation function 
and the critical temperature T^. In particular, the result obtained for u is quite insensitive 
of the estimate of T^, as will be discussed in more detail below. We summarize our esti- 
mates (with error estimates in parentheses): u = 0.80(4), rj = 0.38(3) and = 0.454(4). 
Within the errors, these numbers are in agreement, with a Monte Carlo simulation using 
similar finite-size scaling analysis@ and they are also consistent with other Monte Carlo 
simulation Jii. However, our value of u disagrees with the estimate of = 1 by Thijssen and 
Knops0, who also used the Monte Carlo transfer matrix method, but a different method to 
estimate u; their deviating estimate is likely to be an artifact of their fitting procedure^. 
On the other hand, our value of t] is in very good agreement with their estimate. Moreover, 
our results appear to support the identification of the critical behavior of the fully frustrated 



XY model with the hne of single transitions in XY-Ising model studied recentl; 



II. MONTE CARLO TRANSFER MATRIX 

The Monte Carlo transfer matrix method introduced by Nightingale and Blotelii, is par- 
ticularly useful when, as is the case in the model under consideration, the continuous nature 
of the spin variables does not allow a numerically exact diagonalization of the transfer ma- 
trix. The method is a stochastic version of the well-known power method of calculating the 
dominant eigenvalue of a matrix or integral kernel. More specifically, the dominant eigenvec- 
tor ip^ is approximated by a Monte Carlo time average over weighted walkers representing 
row (or column) configurations. The basic idea is that the (weighted) frequency of occur- 
rence of a particular row (or column) configuration, say S = (si, S2, ■ ■ ■ , Sl), is proportional 
to the magnitude of the corresponding component tpg of the dominant eigenvector. 

The key elements of the algorithm are the following. Walkers are generated in subsequent 
generations labeled by an index t. The generation at time t consists of a sequence of 
a number of rj walkers [(5'i,j, Wi,j), (5*2,4, 1^2,4), (Sr-^^t, tfrt,*)], where the Si^t are row (or 
column) configurations of the form S introduced above, and the Wi^t > are statistical 
weights. This sequence represents a (generally extremely) sparse vector ip with components 

n 

i^s = J2 ^s,s,Wi,u (2) 

i=l 

where 5 is the Kronecker 5-function. One can write 

i^'s' =j:nS'\S)^s = Y.PiS'\S)Dsi^s, (3) 
s s 

where Ds = J2s' T{S'\S) and P(5"|5') = T{S'\S)/ Ds- Since by construction P is a stochastic 
matrix, multiplication by the transfer matrix of a vector of the form can be implemented 
as a stochastic process with transitions from S to 5" with probability P{S'\S). That is, to 
update generation t to t + 1, new walker states are sampled with probability P{St+i\St), 
while in each transition the weight of walker i is updated to Wt^t+i = Dsw^/ct+i- For 
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reasons of efficiency the weights are kept close to unity by duphcating walkers with great 
statistical weights and by eliminating walkers with low weights. To ensure that r^+i remains 
close to its initial value tq, one can choose Ct+i = Xtrt/ro, where is a moving or cumulative 
average estimate of the dominant transfer matrix eigenvalue. Assuming that the generation 
counter t is reset to unity upon equilibration, the largest eigenvalue can be estimated from 
a sequence of T generations as 

\ = j2ct+iWt+jj2Wt, (4) 

t=l t=l 

where Wf = J2iWi^t' denotes the total statistical weight of all walkers of generation t'. It 
should be noted that this estimator of the dominant eigenvalue suffers from a bias due to 
the correlations of the population control constants q with the weights Wt- These correla- 
tions have the effect of suppressing large contributions to the estimator (^) and enhancing 
small oneS. As a result, the estimator tends to underestimate the exact transfer matrix 
eigenvalue, but this effect can be reduced by choosing a large target number of walkers, 
tq. Further algorithmic details, such as how to correct for this bias without increasing the 
average population size and an alternative, sometimes better choice of the population con- 
trol constant q can be found in Refs. |r^, |T^, [T^, ^ In applying this method to the fully 



frustrated XY model, we have performed extensive calculations using, typically, ro = 30, 000 
walkers and 150, 000 Monte Carlo steps which corresponds to 4.5 x 10^ attempts per site. 

To apply this Monte Carlo method it is necessary to sample configurations S' from the 
distribution P{S'\S) for arbitrary S and to evaluate the quantity Ds at each elementary 
Monte Carlo step. One way to do this is to factor the transfer matrix in a way that amounts 
to building up the square lattice by, e.g., ffist adding bonds perpendicular to the transfer 
direction and then adding new sites and bonds along the transfer direction. The algorithm 
given above is then applied to both factors in succession. The ffist of these operations 
corresponds to multiplication by a diagonal matrix. This leaves the row (or column) config- 
urations unchanged while the weight factors are calculated trivially: they consist of one term 
each. The second matrix factor, which adds L new sites, has a direct product form. Sam- 



pling and reweighting problems in this case reduce to simple one-site problems. In principle, 
this method works, but it yields an inefficient Monte Carlo process: the simple factorization 
generates row (or column) configurations with variables that are too weakly correlated in 
the sense that, as it is added to the lattice, each variable is sampled independently from 
a distribution with direct correlations only to one other nearest-neighbor variable. True 
many-site correlations are subsequently recovered via multiplication by strongly fiuctuating 
weight factors. This, however, causes strong fiuctuations in the number of walkers, which 
results in frequent elimination and duplication of walkers and enhances correlations between 
walkers within each generation, which reduces the overall Monte Carlo sampling efficiency. 

The alternative we have used in this paper requires helical boundary conditional. In 
that case the lattice can be constructed by repetition of identical elementary steps, each 
of which adds one site at a time. This implies that the states S' and S in the transition 
matrix P{S'\S) consist of lattice sites all of which coincide with the exception of one site. 
Consequently only a small number of states S' can be reached from any given state S, i.e., 
P is sparse in this case. The degrees of freedom that are added will have direct correlations 
with those at two neighboring lattice sites, one in the horizontal and one in the vertical 
direction for the square lattice. In principle, given an approximate dominant eigenvector, 
one can transform the transfer matrix and incorporate even more correlations from the 
onset, and thus reduce the fiuctuations of the weights and eigenvalue estimate jH^llii. We 
will return to this later on. 

Use of helical boundary yields a more efficient Monte Carlo process, but the method 
has the disadvantage of producing a lattice with a surface that has a step-defect, as is 
unavoidable when one cuts across a screw. One would expect the presence of this defect to 
lead to unnecessary corrections-to-scaling, which may adversely affect finite-size convergence, 
but in practice this appears not to be a serious problem. 

In fact, the method used in this paper requires use of a transfer matrix of somewhat 
more complicated form than the one of Eq. ^ in that this equation has to be replaced by 
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^'^„ = J2T^'\S"\S')T^'\S'\S)^s. (5) 

S',S 

Although this is no real complication -one can simply devise a stochastic process with steps 
alternating in correspondence to the two matrices T'-^^ and T^^-*- it has prevented our use 
of approximate trial eigenvectors to reduce the noise of the stochastic process. In principle, 
this variance reduction scheme works as follows. Eq. ^ can be replaced by the equivalent 
equation 

i^'s' = T.nS'\S)^s (6) 

s 

where T(S"|S') = 75'/T(S"|S')/75' and ip and ip' are similarly similarity transformed. As long as 
the components of 7 are all of one sign and do not vanish, the same Monte Carlo method can 
be applied to T instead of T. This approach satisfies the following zero- variance principle: in 
the ideal limit where 7 is the dominant left eigenvector of T, the largest eigenvalue of T can 
be estimated with vanishing statistical error. More realistically, as this limit is approached 
the variance of the Monte Carlo process decreases. In practical applications, one chooses a 
trial vector 7(^1,^25 ■ ■ ■) which depends on variational parameters Pi. These are optimized 
by minimization of the variance over S of J2s' ls'T{S'\S)/'ys with states selected with prob- 
ability proportional to 7|. This minimization is accomplished approximately by minimizing 
the variance over a relatively small set of states S generated by Monte Carlc@. For a transfer 
matrix of the form of Eq. ^ the same process is still possible in principle, but it is more com- 
plicated. In fact there are two alternatives. One can define a Kronecker-product-like, two-site 
transfer matrix T^^^^S", S'\S', S) = T^^\S"\S')T^^\S'\S) and base the Monte Carlo process 
on the derived stochastic matrix P{S", S'\S', S) oc 'y^g,f]c;iT''^'^\S" , S'\S', S)/'-f^g,'^l, where ide- 
ally 7(12) is the dominant left eigenvector of T^^^\ We note that the pairs {S", S') and (5", S) 
differ at two lattice sites, and therefore this approach requires simultaneous sampling of two 
site variables, which renders the algorithm unnecessarily slow. The alternative is to employ 
-as we have done in this paper- a matrix product of two single-site transfer matrices rather 
than a single two-site transfer matrix, but this approach requires two trial vectors. In terms 
of these, one defines f^^\S'\S) = 7^VT(i)(5'|5)/7f , and f^^]s = jPt'^^^S'IS) /-f^g^ . It is 

7 



straightforward to construct a transformed process which again has a zero variance princi- 
ple. This time it requires that 7^^-' be the dominant left eigenvector of the product matrix 
2"(i)2"{2) that 7*^^-* be the dominant left eigenvector of T^'^^T^^\ Again adjustable param- 
eters of 7*^^) and 7*^^-* can be chosen by minimization of the (appropriately weighted) sum of 
the variances of J2s' ls'T{S'\S)/'ys and J2s' ls'T{S'\S)/'ys- In contrast with the alternative 
of a two-site transfer matrix, the presence of two matrices slows down the algorithm only in 
the initial stage of parameter optimization with this approach. 

III. CHIRAL DOMAIN WALL FREE ENERGY 

For an infinite strip of width L, the reduced free energy / per lattice site can be obtained 
from the largest transfer matrix eigenvalue, X{L,K), via the relation / = — log A, where 
K = J/k^T is the reduced coupling constant. For any given L this quantity, /, depends 
on the choice of boundary conditions. By suitable choice of the latter, as specified in detail 
below, the chiral domain wall free energy can be obtained from the free energy difference, 
denoted by A/. For finite-size analysis, a convenient quantity is the domain wall energy per 
L lattice units of length: 

AF{K,L) = L'AfiK,L). (7) 

Since the fully frustrated XY model is spatially anisotropic, one can devise two different 
types of boundary conditions to compute the domain wall free energy. These two types are 
associated with two different transfer matrices obtained by choosing the transfer direction 
to be either horizontal or vertical, as shown in Fig. |I|. If the transfer direction is horizontal 
(Fig. |l|.a), one is forced to use helical boundary conditions with a pitch of two measured in 
lattice units, so as to match up the vertical antiferromagnetic bonds, which have a periodicity 
of two in the transfer direction. This is the construction used in Ref. 0. As indicated 
in Fig. |I|.a, only strips with L even will match the antiferromagnetic pattern of the local 
chiralities Xp in the ground state; for strips with odd L the boundary conditions will introduce 
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a chiral domain wall along the infinite horizontal direction. Calculation of the domain wall 
energy requires that the boundary conditions be varied at constant L. With the horizontal 
transfer matrix this can be done approximately only, and we chose to use the difference of 
the free energy computed directly for L, and the free energy obtained by linear interpolation 
between sizes L — 1 and L + 1. This still leaves two possibilities, depending on whether L is 
odd or even. As an alternative we also performed calculations using a vertical transfer matrix 
illustrated in Fig. ^ In this case one has a choice between helical boundary conditions with 
a pitch of one or two lattice units. For each value of L precisely one of these two boundary 
conditions will force the presence of a domain wall and this offers a convenient way to 
determine the chiral domain free energy without interpolation. 

Calculations of the free energy using the horizontal and vertical transfer matrices were 
performed as a function of strip width in a very narrow range of couplings AK = 0.04 
around the estimate of the critical coupling obtained by Monte Carlo simulations of Lee et 
a/.0. Note that the range we use is about ten times smaller than the range used in Ref. |10. 
The data for the chiral domain wall free energy using both implementations of the transfer 
matrix are shown in Figs. ^ and ^ 



IV. FINITE-SIZE SCALING ANALYSIS 

To determine the critical temperature and critical exponents, we make use of the following 
finite-size scaling relation for the domain wall free energy 

AF{K, L) = A{L^/^AK), (8) 

where ^4 is a scaling function and AK = K — is the deviation of the coupling constant 
from its critical value Kc. For sufficiently small values of its argument, the scaling function 
can be expanded as 

AF{K, L) = ao + aiL^/^AK + a2{L^''' AKf + • • • , (9) 
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which shows that AF{K, L) is constant as as function of L for K = Kc, but for that value of 
the couphng constant only. This behavior is apparent in Figs. |^ and |^ for at least L > 8 if we 
attribute the deviations for smaller system sizes to corrections to scaling. To the extent that 
the quadratic and higher order terms in Eq. ^ can be neglected, the critical exponent u can 
be obtained from the condition that AF(K, L) be linear in L^^" for fixed AK. In fact, in this 
approximationii, v can be obtained from the slope of a log-log plot of S* = dAF(K, L)/dK 
vs L, which gives l/u . From the data of Figs. ^ and ^, we have obtained S* as a function 
of L for the horizontal and vertical transfer matrix using this procedure. The results are 
indicated in Fig. |[ The slopes of the straight lines in the log- log plot, corresponding to the 
results for vertical and horizontal transfer matrices, agree within the errors, as one would 
expect, providing a check to the consistency of our data. From this plot we can estimate 
u = 0.80(4). It is interesting to note that within the linear approximation -as is the case for 
the free energy barrier in the histogram of chirality in the MC simulation^!, where a similar 
finite size scaling is possible^ - the estimate of the critical temperature = 1 / and the 
estimate of the thermal exponent u are in fact independen'^, but this property is preserved 
only approximately in our generalized finite-size analysis, where we fitted the domain wall 
energies to the form of Eq.(^ with as fitting parameters K^, ly, and with i = 0,1, and 
2, as obtained by truncation of the scaling function beyond second order. Note that the 
correlation function exponent rj can be obtained from the universal amplitude Oq in Eq. ^ 
using the results of conformal invariance in two-dimensionS, = tti]. 

The results of the finite-size scaling based on Eq. ^ fits are summarized in Tab. | and 
the corresponding scaling plots are shown in Figs. || and |[ The table contains estimates of 
the statistical errors associated with the least-squares procedure. In some cases the were 
too large to be attributable to chance and as consequence these statistical error estimates 
are to be treated with suspicion. This is also evident from the discrepancies between the 
various independent estimates. Under the circumstances all we can do is to take the mutual 
differences of various estimates as (admittedly unsatisfactory) error estimates. Thus we 
obtain our estimates: Tc = 0.454(3), u = 0.80(5), and t] = 0.38(2). This result for u 
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is inconsistent with the pure Ising value of z/ = 1, but it is in agreement with the result 
V = 0.85(3) from a similar finite size scaling analysis of the free energy barrier in the 
histogram of chirality from other Monte Carlo simulationJil. Our estimate also agrees 
with a recent estimate u = 0.875(35) based on a finite-size scaling analysis of correlation 
functionJli. We therefore conclude that the estimate of z/ = 1 obtained by Thijssen and 



Recently, the phase diagram of the two dimensional XY-Ising model defined by the 
Hamiltonian 



where a = ±1, has been studied in some detailtllEj. This model is expected to describe the 
critical behavior of a class of systems in which U{1) and Z2 symmetries play a simultaneous 
role, a class of which the fully frustrated XY model is a special case. The critical behavior 
of the square lattice fully frustrated XY model is represented by the behavior of a particular 
(unknown) point in the parameter space {A, C) of this model. The phase diagram of model 
defined in Eq. ^ consists of three branches, in the ferromagnetic region, joining at C ~ 
0. Along one of these branches {C < 0) the transition takes place at a critical point 
with simultaneous criticality of U{1) and Z2 order parameters. Along this line, the critical 
behavior appears to be non- universal with v varying from 0.76 to 0.84 and t] varying from 
0.25 to 0.5 along a segment of the line. Further away from the branch point, this line of single 
transition appears to become first order. Our results for the chiral critical exponents of the 
square lattice fully frustrated XY obtained by the Monte Carlo transfer matrix method do in 
fact appear to be consistent with the critical exponents of a particular point along the line of 
continuous, single transitions of the XY-Ising model. Once this notion is accepted, it implies 
that the fully frustrated XY has a single nonuniversal transition. The non-universality of 
this transition in turn suggests that the critical exponents of the square and triangular lattice 
fully frustrated XY models may differ from each other, although they may be quite closed. 
Besides critical exponents, another important quantity that can be inferred from Monte 




is likely to be an artifact of their fitting procedure. 



(3H 



[A{1 + arCTr') COs(6'r — 6*^') + Ccr^CTr'], 



(10) 



11 



Carlo transfer matrix calculations, is the central charge c, which classifies the possible con- 
formally invariant critical theoriesil. The central charge is related to the amplitude of the 
singular part of the free energy per site, at criticality, in the infinite strip by 

TTC 

f(K,,L) = f^ + — (11) 

which is valid asymptotically for large L. Fitting the data for f{K, L) closest to the estimated 
critical temperature Tc, we obtain c = 1.61(3) from the strips of 8 < L < 14 using both 
horizontal and vertical Monte Carlo data. This result agrees with the estimate of the central 
charge first obtained by Thijssen and KnopJli, also using Monte Carlo transfer matrix 
calculations, and appears to be significantly larger than c = 3/2, which would be expected 
if the transition was single, but decoupledil. However, one cannot be certain of this value 
unless calculations are done at sufficiently large L such that small-L corrections to the above 
asymptotic expression are negligible. Using only the Monte Carlo transfer matrix data for 
6 < L < 14, may not allow us to extrapolate to the L large limit and it is quite likely that 
this estimate of c is subject to systematic errors. In fact, recent large L calculations for the 
related XY-Ising mo deli show a significant decrease of c with increasing L. Unfortunately, 
our data for fully frustrated XY model for L > 14 turn out to be rather noisy. This problem 
can be remedied in principle by using variance reductions techniquesSHS, but as discussed 
this is complicated in the case of the square lattice fully frustrated XY model because of 
the presence of both ferromagnetic and antiferromagnetic bonds. The case of the triangular 
lattice is more straightforward in this respect. 



V. CONCLUSION 

We have studied the finite-size behavior of the chiral domain wall free energy of the 
square lattice fully frustrated XY model on an infinite strip, using the Monte Carlo transfer 
matrix. From a finite-size scaling analysis of data for strip widths up to 14 lattice spacing, 
we have estimated the critical temperature and chiral critical exponents. The latter appear 
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to be significantly different from the pure Ising values and in agreement with other recent 
Monte Carlo simulations. The results are also consistent with the identification of the critical 
behavior of the fully frustrated XY model with that of a point of the line of single transitions 
in XY-Ising model studied recently^. The value of the central charge c = 1.61(3) is found 
to be consistent with the estimate first obtained by Thijssen and Knops, but we cannot rule 
out the possibility that all results obtained for this and similar models are skewed by large, 
slowly decaying corrections-to-scaling. 
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transition due to charging effects in a ladder of Josephson junctions arrays at a half a flux 
quantum per plaquette0. 
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FIGURES 

FIG. 1. Schematic of the process of building up the infinite strip using a transfer matrix. 
This can be implemented using a Monte Carlo transfer matrix along the horizontal direction with 
double helical boundary conditions (a) or along the vertical direction with single (b) or double (c) 
helical boundary conditions. At each step a new configuration is obtained from the previous one, 
indicated by •, by adding a new spin x through a Monte Carlo process and relabeling the sites. 
The lb indicates the antiferromagnetic pattern of the chirality in the ground state. 

FIG. 2. Data for the chiral domain wall free energy /S.F{K, L), obtained from horizontal Monte 
Carlo transfer matrix calculations, in a small range around Tc. 

FIG. 3. Chiral domain wall free energy AF(K, L), obtained from vertical Monte Carlo transfer 
matrix calculations, for the same couplings K of Fig. |2[ 

FIG. 4. S = dAF{K, L) /dK obtained from Monte Carlo transfer matrix data of Figs. |2| and ^. 
The exponent l/v is obtained from the slopes of the linear fit of horizontal (thin line) and vertical 
(thick line) Monte Carlo transfer matrix calculations. 

FIG. 5. Scaling plot of domain wall free energy data obtained for the horizontal transfer matrix: 
AF vs L^^'^AK. The solid curve is the fitted scaling function A, as given by the expansion in Eq. 
The data for L = 6 were not used for the fit. 

FIG. 6. Same as Fig. H for the vertical transfer matrix. 
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TABLES 

TABLE I. Results for critical temperature and critical exponents, obtained from finite-size scal- 
ing analysis of domain wall energies. Standard errors are indicated parenthetically. As explained 
in the text, these standard error reflect only statistical uncertainties, which are presumably con- 
siderably smaller than the errors due to corrections to scaling. The column labeled L indicates 
which system sizes were used in the fits: Li,L2 (AL) stands for sizes from Li to L2 in steps of 
AL. Horizontal and vertical transfer matrix are indicated by v and h under the heading transfer. 
The last column is the per degree of freedom. 



Tc u r) L transfer 

0.4544(3) 0.77(5) 0.378(3) 10, 14 (2) v 11 

0.4555(2) 0.83(4) 0.369(2) 8, 14 (2) v 14 

0.4538(2) 0.83(5) 0.394(1) 8, 12 (2) h 1 
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